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The  rupture  processes  at  a  propagating  crack  tip  and  the  dissipation  of 
mechanical  energy  due  to  inelastic  material  behavior  in  the  highly  deformed 
near- tip  zone,  provide  sources  of  heat  which  can  give  rise  to  substantial 
Increases  in  the  temperature  near  a  crack  tip.  These  temperature  increases 
are  of  obvious  interest  since  they  may  affect  Che  mechanical  behavior  of  the 
near- tip  material,  particularly  its  fracture  properties.  An  early  study  of 
near-tip  temperatures,  motivated  by  interest  in  temperature  effects  on  the 
fracture  toughness  for  fast  fracture,  was  Included  in  a  paper  by  Krafft  and 
Irwin  [1] . 

For  essentially  brittle  materials,  a  propagating  crack  tip  acts  as  a 
sink  of  mechanical  energy.  Parc  of  this  energy  will  be  converted  into 
surface  energy  of  Che  new  crack  faces.  It  is,  however,  generally  assumed 
that  the  remainder  and  perhaps  larger  part,  is  converted  into  heat.  In  the 
mathematical  formulation  a  propagating  crack  edge  is,  therefore,  a 
line  source  of  heat,  and  hence  Che  temperature  at  a  propagating  crack  edge 
will  display  a  logarithmic  singularity,  see  e.g.  Bui  et  al.  [2].  To 
circumvent  this  temperature  singularity,  some  investigators  have  distributed 
the  energy  release  rate  at  the  crack  tip  over  a  small  area.  Studies  of  the 
steady- state  temperature  induced  by  regions  of  near- tip  heat  production 
which  move  with  a  constant  velocity,  have  been  carried  out  by  Weichert  and 
Schdnert  [31-[5].  Other  studies  are  those  of  DOll  [6]  and  Fuller  eC  al  [7], 
Recently  Kuang  and  Atluri  [8]  have  employed  a  finite  element  procedure  to 
calculate  transient  crack-tip  temperatues,  considering  certain  nonlinear 
effects  due  to  temperature  dependent  material  properties  and  convective  and 
radiative  heat  losses . 
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When  Che  crack  propagates  in  an  inelastic  material,  the  dissipation  of 
mechanical  energy  in  the  vicinity  of  Che  crack  tip  provides  another  source 
of  heat.  The  corresponding  transient  crack- tip  temperatures  have  been 
calculated  by  Rice  and  Levy  [9],  for  the  case  that  the  plastic  work  rate  is 
computed  by  Che  use  of  Che  Dugdale  model. 

For  time -dependent  (transient)  temperature  fields,  Che  present  authors 
have  shown  [10]  that  near  the  center  of  a  source  region,  and  under  certain 
restrictions  on  Che  time  of  observation,  Che  effect  of  heat  conduction  on 
the  maximum  temperature  is  negligible  if  a  dimensionless  constant  is 
sufficiently  larger  chan  unity.  The  constant  involves  Che  velocity  and  a 
characteristic  length  parameter  of  Che  source  region,  and  the  thermal 
diffusivity.  When  this  condition  holds,  Che  adiabatic  approximation  which 
completely  ignores  heat  conduction  gives  a  very  acceptable  approximation  to 
Che  temperature  field. 

In  the  present  paper,  Che  adiabatic  approximation  has  been  used  Co 
calculate  the  crack  tip  temperature  for  a  crack  which  propagates  in  a 
viscoplastlc  material.  The  mechanical  behavior  of  the  material  is 
represented  by  a  model  proposed  by  Bodner  and  Partom  [11].  The  general  form 
of  the  near- tip  fields  of  stress  and  deformation  for  this  model  have  been 
studied  by  Achenbach  et  al  [12]  and  Sung  and  Achenbach  [13].  The  Bodner - 
Partom  model  accounts  for  flux  of  energy  into  the  crack  tip,  as  well  as  for 
considerable  dissipation  of  mechanical  energy  in  the  near- tip  region.  In 
this  paper,  Che  temperatures  due  to  the  dissipation,  of  energy  have  been 
analyzed  in  some  detail.  For  a  specified  crack-tip  speed  (300  m/s), 
particular  attention  has  been  devoted  to  the  dependence  of  Che  maximum 
temperature  on  a  characteristic  parameter,  D^,  of  the  Bodner-Partom  model. 
This  parameter  defines  the  asymptotic  limit  of  Che  inelastic  strain  rate  as 
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the  stress  increases.  The  maximuni  temperature  Increases  with  increasing 

values  of  D. . 

o 

2.  Constitutive  Model 

In  this  paper  we  use  the  constitutive  relations  for  a  viscoplastic 
material  proposed  by  Bodner  and  Partom  [11].  The  basic  equation  that 
distinguishes  the  Bodner-Partom  model  from  other  constitutive  models,  is  the 
one  which  defines  the  rate  of  inelastic  (plastic)  strain  as: 

‘ij  "  °o  2  .  (2.1) 

CJ  2^ 


where  s^j  are  the  components  of  thu  stress  deviator, 


®ij  “  ‘^ij  *  3  ^ij'^kk  ' 


(2.2) 


and 


•^2  "  2  ®ij®ij 


(2.3) 


Also 


-mW 

Z  -  +  (Z^  -  Z^)e  P  (2.4) 

In  Eqs.  (2.1)  and  (2.4),  Z^,Zi ,m,n  and  D  are  material  constants,  and  W  is 

O  1  o  p 

the  work  over  the  inelastic  strains 

(2.5) 


3 


The  Rodner>Partom  constitutive  equations  have  the  convenient  property 
that  no  separate  specification  of  a  yield  criterion  Is  required,  and  loading 
and  xinloadlng  are  governed  by  the  same  set  of  equations.  Both  elastic  and 
Inelastic  deformations  are  present  at  all  stages  of  loading  and  unloading. 
For  loading  In  uniaxial  tension  It  Is  simple  to  plot  the  Inelastic  strain 
rate,  versus  Figure  1  shows  that  remains  very  small  at 

small  stresses,  with  a  rapid  change  over  a  small  range  of  (the  "yield" 

region) ,  whereupon  an  asymptotic  limit  (proportional  to  0^)  Is  approached  as 
becomes  large. 

It  has  been  shown  by  Sung  and  Achenbach  [13]  that  a  bounded  Inelastic 
strain  rate  Implies  a  square-root  singular  near- tip  stress  field.  In  the 
moving  coordinates  we  have; 


‘'ll 

‘'IZ 

"22 


K^(t) 


3  115 

J  cosj^  +  *  cos2^ 

-  ^  sinj^  +  ^  sin2^ 

^  cosi^tf  -  ^  cos^^ 

4  2  4  2 


(2.6) 


It  should  be  noted  that  the  usual  factor  (2ir)  In  the  denominator  has  been 
omitted,  and  that  the  variations  with  angle  are  the  usual  ones  for  elastic 
near- tip  stress  fields. 

Sung  and  Achenbach  [13]  have  used  the  finite  element  method  to 
calculate  the  quasi-static  stress  Intensity  factor,  K^(t).  They  considered 
a  case  of  plane  strain  and  small-scale  yielding,  for  a  body  containing  a 
stationary  crack  which  Is  loaded  at  time  t  -  0.  It  was  assumed  that  the 
load  Immediately  generates  the  Initial  elastic  field  everywhere  In  the 
cracked  body.  In  particular  the  Initial  fields  near  the  crack  tip  were 
assumed  to  consist  of  the  usual  static  elastic  singular  fields.  At  time 
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t  •  0^  the  crack  starts,  however,  to  propagate  with  a  constant  velocity. 

^  £  E 

Curves  for  Kj(t  )/Kj,  where  Kj  Is  the  initial  elastic  stress  intensity 
factor,  are  shown  in  Fig.  2,  for  various  values  of  the  crack  tip  speed.  The 
dimensionless  time  is  defined  by 


t*  -  t(/iD^)/(Z/y3) 


(2.7) 


For  the  purpose  of  the  present  paper  the  curves  of  Fig.  2  are  fitted  by 


(2.8) 


where  a  and  a  depend  on  the  crack- tip  speed. 

It  should  be  noted  that  the  components  of  ,  as  defined  by  Eq.(2.1), 

are  bounded  at  a  crack  tip.  Even  though  s^j  blows  up  as  a  crack  tip  is 

approached,  in  the  term  Sj^j/(J2)  that  singularity  is  cancelled  by  the  one 
H  n 

of  (Jj)  .  The  term  exp[ - (Z*/3J2)  /2]  is  bounded  whether  J2  is  singular  or 
not.  At  a  crack  tip  we  therefore  have 


JP  - 
ij 


D  s.  . 
9 


^2®ki®ki^ 


(2.9) 


This  expression  displays  the  strong  dependence  of  the  near- tip  inelastic 
strain  rate  on  the  parameter  D^. 

Eqtiations  (2.5),  (2.6)  and  (2.8)  are  now  used  to  compute  the  rate  of 
dissipation  of  mechanical  energy  in  the  immediate  vicinity  of  the  crack  tip. 
The  result  is 
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(2.10) 


*p(x,y,t) 


,  „E  -at  TICS') 


where 

S(tf)  -  2{^[1  +  4v(i/-l)]cos*^tf  +  J  (2.11) 

3.  Adiabatic  Approximation 

The  calculation  of  the  temperature  at  a  propagating  crack  tip  is  a 
moving  heat-source  problem.  The  rupture  processes  right  at  the  propagating 
tip  and  the  dissipation  of  mechanical  energy  in  the  vicinity  of  the  crack 
tip,  provide  the  moving  heat  sources. 

When  the  source  is  distributed  and  of  the  general  form  Q^f(X,Y,t), 
where  X  and  Y  are  the  moving  coordinates  shown  in  Fig.  1,  and  if  the  system 
of  sources  moves  at  a  constant  velocity,  V,  the  temperature  in  the  moving 
coordinate  system  may  be  written  as,  see  (10], 


T(X,Y,t)  -  ,  _  /  f(?,n.s)  exp(- 

°  t-s  A 


( .(X- f ItV-I.?: -  s }J  -LY-jH) 

2(t-i) 


(3.1) 


where  the  dimensionless  variables  are  defined  as 


(X,Y)  -  (X,Y)/i  ,  (i.n)  -  .  (t,s)  -  (2«/f2)(t,s),  (3.2a,b,c) 


T  -  T/(Q^f*/8k)  ,  V  -  (f/2»c)V. 


(3.3a,b) 
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In  these  expressions  i  is  a  characteristic  length  of  the  area  A,  k  Is  the 
thermal  conductivity,  and  k  is  the  thermal  diffusivlty 


*  -  k/pc  ,  (3.4) 

where  p  is  the  mass  density  and  c  is  the  specific  heat. 

Equation  (3.1)  reveals  the  potentially  dominant  influence  of  V. 

Suppose  that  X.Y  <  1,  and  that  V  »  1.  The  exponential  term  will  then  be 

very  small,  unless  t-s  «  1.  When  t-s  «  1  the  terms  [ (X-^)+V( t- s) ] *72 ( t- s) 

and  (Y-f;)2/2(t-s)  will,  however,  be  very  large,  and  the  exponential  term 

will  decay  rapidly  as  ?  and  t)  vary  from  X  and  Y.  These  observations  suggest 

that  for  large  values  of  V,  the  major  contribution  from  f(|^,^,s)  comes  from 

the  neighborhood  of  ?-X+V(t-s) ,»/-Y.  Hence  we  bring  f(X,Y,s)  outside  the 
integration  over  A.  Since  the  exponential  term  is  either  very  small 

(t-s  -  0(1))  or  decays  rapidly  (t-s  «  1),  the  limits  of  integration  may  be 
extended  to  ±<o.  Use  of  the  integral 


then  finally  reduces  Eq.(3.1)  to 


T(X,Y,t)  -  4  f(X+V(t-I),Y,s)ds  (3.6) 

o 
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It  follows  from  (3.3b)  that  V  is  large  for  Vi  large  and/or  k  small. 
However,  even  for  moderate  values  of  Vi,  one  may  find  quite  large  values  of 

V.  For  example,  for  steel  we  have 

pc  -  3.64  X  10«  J/m»  <»C,  k  -  5.36  X  10  J/ms  “C  (3.7a,b) 

and  hence 

K  -  ^  -  1.47  X  10'^  mS/s  (3.8) 

pc.  ' 

It  then  follows  from  Eqs.(3.3b)  and  (3.8)  that 

V  -  Vi  X  lOS/2.94  (3.9) 


Clearly  V  can  be  quite  large  even  for  small  values  of  V  or  i . 

For  certain  simple  moving  source  regions  the  temperatures  according  to 

Eqs.  (3.1)  and  (3.6)  have  been  compared  in  Ref. [10].  For  V  -  50  the 
agreement  is  very  good  at  small  times ,  but  it  degenerates  as  time  increases . 

Equation  (3.6)  may  be  called  an  adiabatic  approximation.  It  assumes 
that  the  source  region  moves  sufficiently  fast  that  the  effect  of  heat 
conduction  is  negligible  in  the  source  region,  at  least  at  relatively  small 
times.  The  purely  adiabatic  approximation  can  also  be  obtained  directly  by 
simply  setting  k  0  in  the  heat  conduction  equation.  The  equation  for  the 
temperature  then  reduces  to 

pet  -  Q  (3.10) 
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Suppose  Q  -  QQf(x  -  a(t),y,t)  ts  given  in  the  stationary  coordinate  system. 
The  temperature  at  a  fixed  point  x,y  relative  to  the  stationary  coordinate 
system  is 

Q  t 

T(x,y,t)  -  -7  /  f(x-a(s) ,y,s)ds  (3.11) 

^  o 

where  we  assume  T(x,y,t)  -  0  for  t  <  0.  Equation  (3.11)  is  generally  quite 
easy  to  evaluate,  even  numerically  when  f(x  -  a(t),y,t)  is  known  only  in  the 
form  of  numerical  data.  In  the  moving  coordinates  Eq.  (3.11)  becomes 

Q  t 

T(X,Y,t)--7  /  f(X+a(t)-a(s),Y.s)ds  (3.12) 

pc  J 

The  corresponding  dimensionless  results  in  the  moving  coordinate  system  are 
obtained  by  using  the  dimensionless  quantities  defined  by  Eq.(3.2)  and 
Eq.(3.3).  Note  that  the  actual  value  of  k  appears  in  these  quantities,  even 
though  Eq.(3.12)  is  for  k  -  0.  We  have 


t  _  _ _ _ 

T(X,Y,t)  -4/  f(X+a(t)-a(s) ,Y,s)ds  (3.13) 

o 

In  the  absence  of  heat  conduction,  and  since  the  source  region  moves 
parallel  to  the  x  axis,  the  coordinate  y  simply  appears  as  a  parameter.  It 
is  noted  that  for  a  source  moving  with  constant  velocity  V,  Eq.(3.13)  is 
identical  to  Eq.(3.6). 


f 


rr 


»  i*.  'M 


Under  the  assumption  that  the  dissipated  mechanical  energy  is  converted 

into  heat,  the  temperature  generated  by  the  singular  part  o£  Up(X,Y,t)  can 
now  be  obtained  by  substitution  of  (2.10)  into  Eq.(3.12}.  The  result  is 


T(X,Y.t) 


D 

o  I  _ 

pc 


t 

; 

o 


-as  . 

e  -d£ 


([X+a(t)-a(s)]2  +  y2)^/^ 


(3.14) 


where  a  is  defined  by  Eq.(2.8).  For  constant  crack  tip  speed,  V,  and  after 
some  manipulation,  we  find  on  the  crack  line  (Y  -  0) ; 


T(X,0,t) 


D^K^(O) 

pc 


^-a(t+X/V) 

(aV)^  g(aX/V,ot) 


(3.15) 


where 


a(X/V+t)  s 

g(aX/V,eit)  -  /  ~Tn 

oX/V  s' 


At  the  crack  tip  the  temperature  further  simplifies  to 


(3.16) 


T(0,0,t) 


DqK^(O) 

pc 


-at 


e 


g(0,ot) 


(3.17) 


The  function  e  g(0,Qt),  which  is  known  as  Dawson's  integral,  has  been 
plotted  in  Fig.  3.  It  is  noted  that  the  maximum  value  of  the  integral  is 
approximately  1.1.  Hence,  at  the  crack  tip  we  have 


T 

max 


-Li 

(«V)‘* 


D  K^(0) 
g  I 
pc 


(3.18) 
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It  Is  of  Interest  to  plot  the  dependence  of  T  on  D  .  To  do  that,  we 

mAX  o 

use  the  nondimens lonal  time  introduced  by  Eq.(2.7).  The  crack- tip  speed,  V, 
is  then  related  to  the  normalized  crack-tip  speed,  V*,  by  the  following 
expression: 


V 


z» 


373  V* 


(3.19) 


Let  us  consider  the  following  numerical  values 
Kj  -  70  MPa  m'*  -  0.7  X  10«N/m^/^ 

H  —  .44  X  10®  N/mm^  ,  —  0.2 


Z  -  1550  N/mm2  ,  V  -  300  m/s 


Equation  (3.19)  then  becomes 

D  V*  -  10'^  1/s  . 

o  ' 


(3.20) 


Therefore,  for  various  values  of  D^,  we  can  compute  the  corresponding  values 
for  V  by  the  above  formula.  The  results  are  shown  in  Table  1.  As  pointed 

'if 

out  earlier,  the  quasi-static  stress  intensity  factor  Kj(t  )  computed 
numerically  in  Ref. [13]  for  various  normalized  crack-tip  speeds  V  (see  Fig. 
2),  can  be  fitted  for  short  time  by  Eq.(2.8),  by  adjusting  the  parameter  a 

'if  'if  ' 

(or  a  if  the  normalized  time  scale  t  -  /iD^/Z/73  is  used) ,  where 
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(3.21) 


it 

For  various  normalized  crack* tip  speeds,  V  ,  the  fitted  values  for  the 
parameter  a  are  shovm  in  Table  1. 

For  a  specific  value  of  V  (in  our  case  V  —  300  m/s)  the  calculation  of 

T  as  a  function  of  D  now  proceeds  as  follows.  For  various  selected 
max  o 

values  of  D^,  the  corresponding  V  are  obtained  from  (3.20).  Next,  the 

corresponding  values  of  a*  are  obtained  by  interpolation  from  Table  1.  For 

each  the  value  of  a  then  follows  from  Eq.(3.21).  Finally  for  each  the 

corresponding  a  is  substituted  in  Eq.(3.17),  to  yield  The  curve  of 

T  versus  D  is  shown  in  Fig.  4.  It  is  noted  that  T  increases 
max  o  °  max 

monotonically  with  D^.  The  increase  is  more  dramatic  at  higher  values  of 
D^.  For  example,  for  -  10®  1/s  we  have  T^^^  »  360°C,  while  for 
D  -  10®  1/s,  we  have  T  =«  20®C. 

Crack  tip  temperatures  due  to  flux  of  energy  into  the  crack  tip  have 
been  computed  in  Ref. [10], 
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Time  varying  stress-intensity  factor,  k(t  )  ■  K^(t  J/K^, 
versus  the  dimensionless  time,  for  various  crack-tip  speeds, 
where  V*  Is  defined  by  Eq.(3.19). 
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